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We present approximation algorithms for almost all variants of the multi-criteria 
traveling salesman problem (TSP), whose performances are independent of the 
number k of criteria and come close to the approximation ratios obtained for TSP 
with a single objective function. 

We present randomized approximation algorithms for multi-criteria maximum 
traveling salesman problems (Max-TSP). For multi-criteria Max-STSP, where the 
edge weights have to be symmetric, we devise an algorithm that achieves an approx- 
imation ratio of 2/3 — e. For multi-criteria Max-ATSP, where the edge weights may 
be asymmetric, we present an algorithm with an approximation ratio of 1/2 — e. 
Our algorithms work for any fixed number k of objectives. This improves over the 
previously known algorithms, which achieve ratios oil/k — e for multi-criteria Max- 
STSP and l/{k + 1) — e for multi-criteria Max-ATSP, and is close to the currently 
best approximation ratios for Max-STSP and Max-ATSP with a single objective 
(61/81 3/4 and 2/3, respectively). To get the approximation ratios of 2/3 — e 
and 1/2 — e, respectively, we introduce a decomposition technique for cycle covers. 
These decompositions are optimal in the sense that no decomposition can always 
yield more than a fraction of 2/3 and 1/2, respectively, of the weight of a cycle 
cover. Furthermore, we present a deterministic algorithm for bi-criteria Max-STSP 
that achieves an approximation ratio of 61/243 ~ 1/4. 

Finally, we present a randomized approximation algorithm for the asymmetric 
multi-criteria minimum TSP with triangle inequality (Min-ATSP). This algorithm 
achieves a ratio of logn + e. For this variant of multi-criteria TSP, this is the 
first approximation algorithm we are aware of. If the distances fulfil the 7-triangle 
inequality, its ratio is 1/(1 — 7) -|- e. 



1 Multi-Criteria Traveling Salesman Problem 
1.1 Traveling Salesman Problem 

The traveling salesman problem (TSP) is one of the most famous combinatorial optimization 
problems. Given a graph, the goal is to find a Hamiltonian cycle of maximum or minimum 
weight (Max-TSP or Min-TSP). 

*To be presented at the 26th Int. Symposium on Theoretical Aspects of Computer Science (STAGS 2009). 
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An instance of Max-TSP is a complete graph G = iy,E) with edge weights w : E ^ 
The goal is to find a Hamiltonian cycle of maximum weight. The weight of a Hamiltonian cycle 
(or, more general, of any set of edges) is the sum of the weights of its edges. If G is undirected, 
we have Max-STSP (symmetric TSP). If G is directed, we obtain Max-ATSP (asymmetric 
TSP). 

An instance of Min-TSP is also a complete graph G with edge weights w that fulfil the 
triangle inequality: w{u,v) < w{u,x) + w{x,v) for all u,v,x G V. The goal is to find a 
Hamiltonian cycle of minimum weight. We have Min-STSP if G is undirected and Min-ATSP 
if G is directed. In this paper, we only consider the latter. If we restrict the instances to 
fulfil the 7-triangle inequality {w{u,v) < 7 • {w{u,x) +w{x,v)) for all distinct u,v,x &V and 
7 G [i, 1)), then we obtain Min-7-ATSP. 

All variants introduced are NP-hard and APX-hard (Min-7-ATSP is hard for 7 > ^). Thus, 
we have to content ourselves with approximate solutions. The currently best approximation 
algorithm for Max-STSP achieves an approximation ratio of 61/81 [10], and the currently best 
algorithm for Max-ATSP achieves a ratio of 2/3 [23] . Min-ATSP can be approximated with a 
factor of I • log2 n, where n is the number of vertices of the instance [14]. Min-7-ATSP allows 
for an approximation ratio of minj^^, 2 J""*^ 3 } [7,9]. 

Cycle covers are often used for designing approximation algorithms for the TSP [4, 5, 7, 9, 
10, 14,23,24]. A cycle cover of a graph is a set of vertex-disjoint cycles such that every vertex 
is part of exactly one cycle. The general idea is to compute an initial cycle cover and then to 
join the cycles to obtain a Hamiltonian tour. This technique is called subtour patching [17]. 
Hamiltonian cycles are special cases of cycle covers that consist of a single cycle. Thus, the 
weight of a maximum- weight cycle cover bounds the weight of a maximum- weight Hamiltonian 
cycle from above, and the weight of a minimum-weight cycle cover is a lower bound for the 
weight of a minimum-weight Hamiltonian cycle. Moreover, in contrast to Hamiltonian cycles, 
cycle covers of maximum or minimum weight can be computed efficiently by reduction to 
matching problems [1]. 

1.2 Multi-Criteria Optimization 

In many optimization problems, there is more than one objective function. This is also the case 
for the TSP: We might want to minimize travel time, expenses, number of flight changes, etc., 
while maximizing, e.g., the number of sights along the way. This leads to /c-critcria variants 
of the TSP (fe-C-Max-STSP, fc-C-Max-ATSP, fc-C-Min-STSP, fe-C-Min-ATSP for short; if the 
number of criteria does not matter, we will also speak of MC-Max-STSP etc.). 

With respect to a single criterion the term "optimal solution" is well-defined. However, if 
several criteria are involved, there is no natural notion of a best choice. Instead, we have to 
be content with trade-off solutions. The goal of multi-criteria optimization is to cope with this 
problem. To transfer the concept of optimal solutions to multi-criteria optimization problems, 
the notion of Pareto curves (also known as Pareto sets or efficient sets) was introduced (cf. 
Ehrgott [12]). A Pareto curve is a set of solutions that can be considered optimal. 

We introduce the following terms only for maximization problems. After that, we briefly 
state the differences for minimization problems. 

An instance of fc-C-Max-TSP is a complete graph G with edge weights wi, . . . ,Wk : E ^ Q+. 
A Hamiltonian cycle H dominates another cycle H' if Wi{H) > Wi{H') for all i E [k] = 
{1,. . . ,k} and Wi{H) > Wi{H') for at least one i. This means that H is strictly preferable 
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io H'. A Pareto curve of solutions contains all solutions that are not dominated by another 
solution. For other maximization problems, /c-criteria variants are defined analogously. 

Unfortunately, Pareto curves cannot be computed efficiently in many cases: First, they 
are often of exponential size. Second, because of straightforward reductions from knapsack 
problems, they are NP-hard to compute even for otherwise easy optimization problems. Third, 
the TSP is NP-hard already with only one objective function, and optimization problems do 
not become easier with more objectives involved. Therefore, we have to be satisfied with 
approximate Pareto curves. 

For simpler notation, let w{H) = (^wi{H), . . . ,Wk{H)Y Inequalities are meant component- 
wise. A set V of Hamiltonian cycles of V is called an a approximate Pareto curve for (G, w) 
if the following holds: For every Hamiltonian cycle H' , there exists a Hamiltonian cycle H G 
V with w{H) > aw{H'). We have a < 1, and a 1 approximate Pareto curve is a Pareto 
curve. (This is not precisely true if there are several solutions whose objective values agree. 
Furthermore, with this definition, the set of all feasible solutions forms a Pareto curve. If 
one allows succinct representations of sets, then many NP-hard multi-criteria problems are 
trivially solvable in polynomial time. However, our algorithms output sets in a "natural" 
representation, which allows us to remove dominated solutions from the sets easily. Thus, in 
our case this subtlety is inconsequential, and we will not elaborate on it for the sake of clarity.) 

An algorithm is called an a approximation algorithm if, given G and w, it computes an a 
approximate Pareto curve. It is called a randomized a approximation if its success probability 
is at least 1 /2. This success probability can be amplified to 1 — 2""* by executing the algorithm 
m times and taking the union of all sets of solutions. (We can also remove solutions from this 
union that are dominated by other solutions in the union, but this is not required by the 
definition of an approximate Pareto curve.) Again, the concepts can easily be transfered to 
other maximization problems. 

Papadimitriou and Yannakakis [26] showed that (1 — e) approximate Pareto curves of size 
polynomial in the instance size and 1/e exist. The technical requirement for the existence 
is that the objective values of all solutions for an instance X are bounded from above by 
2pW for some polynomial p, where is the size of X. This is fulfilled in most optimization 
problems and in particular in our case. However, they only prove the existence, and for many 
optimization problems it is unclear how to actually find an approximate Pareto curve. 

A fully polynomial time approximation scheme (FPTAS) for a multi-criteria optimization 
problem computes (1 — e) approximate Pareto curves in time polynomial in the size of the 
instance and 1/e for all e > 0. Papadimitriou and Yannakakis [26] showed that multi-criteria 
minimum- weight matching admits a randomized FPTAS, i. e., the algorithm succeeds in com- 
puting a (1 — e) approximate Pareto curve with constant probability. This randomized FPTAS 
yields also a randomized FPTAS for the multi-criteria maximum-weight cycle cover prob- 
lem [25]. 

To define Pareto curves and approximate Pareto curves also for minimization problems, in 
particular for MC-Min-STSP and MC-Min-ATSP, we have to replace all ">" and ">" above 
by "<" and "<". Furthermore, a approximate Pareto curves are then defined for a > 1 and 
an FPTAS has to achieve an approximation ratio of 1 -|- e. There also exists a randomized 
FPTAS for the multi-criteria minimum-weight cycle cover problem. 
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1.3 Related Work 



A considerable amount of research has been done on multi-criteria TSP. Most work is about 
heuristics for finding approximate solutions without any worst-case guarantee, e.g., simulated 
annealing [18], tabu search [20], or genetic algorithms [22]. There are also algorithms for finding 
the exact Pareto curve, although this takes exponential time [8,28]. Gupta and Warburton [19] 
used local search and Langrangean methods to find locally optimal solutions for the so-called 
Tchebycheff approach, where the aim is to find a solution that minimizes maxjgjj^.] • Wi{H). 
Paquete and Stiitzle [27] used local search and Fischer and Richter [15] used dynamic pro- 
gramming methods to solve bi-criteria TSP. 

Angel et al. [2, 3] considered Min-STSP restricted to edge weights 1 and 2. They analyze a 
local search heuristic and prove that it achieves an approximation ratio of 3/2 for A; = 2 and 
of ^ for > 3. 

Ehrgott [11] considered a variant of MC-Min-STSP, where all objectives are encoded into a 
single objective by using some norm. He proved approximation ratios between 3/2 and 2 for 
this problem, where the ratio depends on the norm used. 

Manthey and Ram [25] designed a (2 -|- e) approximation algorithm for MC-Min-STSP and 
an approximation algorithm for MC-Min-7-ATSP, which achieves a constant ratio but works 
only for 7 < l/v/3 ^ 0.58. They left open the existence of approximation algorithms for 
MC-Max-STSP, MC-Max-ATSP, and MC-Min-ATSP. 

Blaser et al. [6] devised the first randomized approximation algorithms for MC-Max-STSP 
and MC-Max-ATSP. Their algorithms achieve ratios oi\+£ for /c-C-Max-STSP and ^ + e 
for fc-C-Max-ATSP. They argue that with their approach, only approximation ratios of j^j^Q^i^^ 
can be achieved. Nevertheless, they conjectured that approximation ratios of Q.{l/\ogk) are 
possible. 

For an overview of the literature about multi-criteria optimization, including multi-criteria 
TSP, we refer to Ehrgott and Gandibleux [13]. 



1.4 New Results 

We devise approximation algorithms for MC-Max-STSP, MC-Max-ATSP, and MC-Min-ATSP. 
The approximation ratios achieved by our algorithms are independent of the number k of 
criteria, and they come close to the best approximation ratios known for Max-STSP, Max- 
ATSP, and Min-ATSP with only a single objective function. Our algorithms work for any 
number k of criteria. 

First, we solve the conjecture of Blaser et al. [6] affirmatively. We even prove a stronger 
result since the performance ratios of our algorithms are independent of k: For MC-Max-STSP, 
we achieve a ratio of 2/3 — e, while for MC-Max-ATSP, we achieve a ratio of 1/2 — e (Section 4). 
Already for k = 2, this is an improvement from ^ — e to | — e for 2-C-Max-STSP and from \ — £ 
io\-e for 2-C-Max-ATSP. The general idea of our algorithm is sketched in Section 2. After 
that, we introduce a decomposition technique in Section 3 that will lead to our algorithms. 
The running-time of our algorithms is polynomial in the input size for any fixed e > and any 
fixed number k of criteria. 

Furthermore, as a first step towards deterministic approximation algorithms for MC-Max- 
TSP, we devise an approximation algorithm for 2-C-Max-STSP that achieves an approximation 
ratio of 61/243 > 1/4. As a side effect, this proves that for 2-C-Max-STSP, there always exists 
a single Hamiltonian cycle that already is a 1/3 approximate Pareto curve. For completeness. 
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we show that this does not hold for fc-C-Max-STSP for k > S, for MC-Max-ATSP, or for 
MC-Min-TSP. 

Finally, we devise the first approximation algorithm for MC-Min-ATSP (Section 6). In 
addition, our algorithm improves on the algorithm for MC-Min-7-ATSP by Manthey and 
Ram [25] for 7 > 0.55 (see Figure 4 on page 24), and it is the first approximation algorithm 
for MC-Min-7-ATSP for 7 G [-^, l) . The approximation ratio of our algorithm is log n + e for 

MC-Min-ATSP, where n is the number of vertices. Furthermore, it is a + e approximation 
for MC-Min-7-ATSP for 7 G [^jl)- Our algorithm is randomized and its running-time is 
polynomial in the input size and in 1/e for any fixed number of criteria. 

2 Outline and Idea for MC-Max-TSP 

For Max-ATSP, we can easily get a 1/2 approximation: We compute a maximum- weight cycle 
cover, and remove the lightest edge of each cycle. In this way, we obtain a collection of paths. 
Then we add edges to connect the paths, which yields a Hamiltonian cycle. For Max-STSP, 
this approach yields a ratio of 2/3 since the length of every cycle is at least three. 

Unfortunately, this does not generalize to multi-criteria Max-TSP, even though {1 + e) ap- 
proximate Pareto curves of cycle covers can be computed in polynomial time. The reason is 
that the term "lightest edge" is usually not well defined: An edge that has little weight with 
respect to one objective might have a huge weight with respect to another objective. Based on 
this observation, the basic idea behind our algorithms is the following case distinction: First, 
if every edge of a cycle cover is a light-weight edge, i.e., it contributes only little to the overall 
weight, then removing one edge docs not decrease the total weight by too much. We can choose 
the edges to be removed such that no objective loses too much of its weight. 

Second, if there is one edge that is very heavy with respect to one objective (a heavy-weight 
edge), then we take only this edge from the cycle cover. In this way, we have enough weight 
for one objective, and we proceed recursively on the remaining graph with k — 1 objectives. 

In Section 3, wc deal with the first case. This includes the definition of when we call an 
edge a light-weight edge. In Section 4, we present our algorithm, which includes the recursion 
in case of a heavy-weight edge. The approximation ratios that we achieve come close, i.e., up 
to an arbitrarily small additive £ > 0, to the 1/2 and 2/3 mentioned above for mono-criterion 
Max-ATSP and Max-STSP. 

3 Decompositions 

3.1 Asymptotically Optimal Decompositions 

Prom any collection P of paths, we obtain a Hamiltonian cycle just by connecting the endpoints 
of the paths appropriately. Assume that we are given a cycle cover C. If we can find a collection 
of paths P (1 C (by removing one edge of every cycle of C) with w(P) > a ■ w(C) for some 
a G (0, 1], then this would yield an approximate solution for Max-TSP. 

Let us define decompositions more formally. Let a G (0,1], and let C be a cycle cover. 
Blaser ct al. [6] called a collection P C C of vertex-disjoint paths an a-decomposition of C if 
w{P) > aw{C). 

Not every cycle cover possesses an a-decomposition for every a. Let k > 1 he the number 
of criteria. Blaser et al. defined G [0, 1] to be the maximum number such that the following 
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holds: every directed cycle cover C with edge weights w = {wi, . . . ,Wk) that satisfies w{e) < 
af ■ w{C) for all e G C possesses an a^-decomposition. The value € [0, 1] is analogously 
defined for undirected cycle covers. 

We have af = ^ and a" = |. We also have > and < a^_i as well as < af._i. 

Blaser ct al. [6] proved af. > and > p Furthermore, they proved the existence 
of log A;)-dccompositions, i.e., 0(f,a^ G Q(l/\ogk), which led to their conjecture that 
ri(l/log/c) approximation algorithms might exist. However, their approximation algorithms 
do not make use of the $1(1/ log fc) decompositions, and they only achieve ratios of ^ — £ for 
A;-C-Max-STSP and — e for A;-C-Max-ATSP. In fact, they indicate that approximation 
ratios of i^_^_o(i) the best that can be proved using their approach. 

For completeness, we make their decomposition result more precise with the next theorem. 
In particular, we show that Q;^,a^ G 0(1/ log /c), which proves that better approximations 
require a different decomposition technique. The new decompositions will be introduced in 
Section 3.2. 

Theorem 3.1. For all 1 < k eN, we have 

1 1 „ 1 1 , 

< or < -r, r-. r ~ — ; ; t ana 



0.78 ■ log2 A; + §~|-lnfc + |- ^ [log^ k\+l " 0.63 • logg k + 1 

« < < . 

1.39-log2A; + 4 2-lnfe + 4- " Llog2 + 2 

Proof. We consider the upper bounds first. It suffices to prove af < i^ilg^ ^ for k = 2^ and 

q" < i^ilg fc for k = 3^ for ^ G N since > a'^_^^ and > ctf+i- Let us first analyze af. for 

k = 2^ For ^ = 0, 1, we already know af < 1/2 = 2+155-3 < 1/3 = 2+^^ [6]. So let 

£>2. 

We now describe a cycle cover C with w{e) < 1 for all e that does not possess a decomposition 
P with w{P) > 1. This proves < l/wi{C) for i G [k]. 

Let e > be a small number that we will specify later on. We hold a "knockout tournament" 
for the k = 2^ objectives. We model a match of / C [k] versus J C [k] by a cycle (e, /) with 
Wi(e) = 1 — £ foT i e I, Wi{e) = for i ^ /, Wi{f) = 1 — e for z G J, and Wi{f) = for 
i ^ J. Winning a match thus yields weight 1 — e as a prize. In the first round, {1} plays 
against {2}, {3} plays {4}, . . . , and {k ~ 1} plays {k}. In the second round, {1, 2} plays {3, 4}, 
{5, 6} plays {7, 8}, and so on. In the "semi-final", which is the (i — l)th round, {1, . . . , A:/4} 
plays {k/A + 1, . . . , k/2} while {k/2 + 1, . . . , 3A;/4} plays {3A;/4 + l,...,k}. The final are three 
matches were {1, . . . , A;/2} plays against {k/2 + 1, . . . ,k}. A complete tournament for ^ = 3 is 
shown in Figure 1(a). 

We call an objective i a p-underdog if i has lost all its matches up to the pth round. 

Lemma 3.2. For p G [f] and for every j G [2^~^], there is a p-underdog in {{j — 1) • 2^' + 1, (j — 
l)-2P + 2,...,j-2P}. 

Proof. The proof is by induction on p. For p = I, this is true since {(j — 1) • 2 + 1} = {2j — 1} 
plays {2j} and only one of them can win. For p > 1, we have one {p — l)-underdog in all 
{{j — 1) ■ 2P~^ + 1,... ,j ■ 2P~^} for j G [2^~^+^]. In every match of the pth round, a team 
{{j - 1) • 2P-^ + 1, . . . , J • 2P-1} plays a team {j ■ 2^'^ + 1, . . . , (j + 1) • 2^"^} for an odd j. Either 
team contains a (p — l)-underdog. One of the (p — l)-underdogs loses again and becomes a 
p-underdog. □ 
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(a) 






1 4 



(b) 



1 4 



Figure 1: (a) Directed tournament for 8 = 2^ objectives. The labels represent the non-zero 
entries of the corresponding vectors. For instance, label 3, 4 represents the vector 
(0, 0, 1 — £, 1 — £, 0, 0, 0, 0). (b) Undirected tournament for 9 = 3^ objectives. 



Thus, before the three final matches, wc have two (i — l)-underdogs. In the now following 
three final matches, these two (i — l)-underdogs play each other. Thus, only one of them can 
achieve the necessary two victories. 

The tournament results stand in one-to-one correspondence to decompositions of the cycle 
cover C that represents it: We put all "winning" edges into P. We conclude that for every 
decomposition P of C, there is an i with Wi{P) < 1 — e < 1. Every objective is involved into 



i + 2 matches. Thus, Wi{C) = {l + 2)-{l-e) and Wi{e 
does not posses a (£^2)\i-£) -decomposition. Hence "'^ 



< 1 



(^+2).(l-e)' 



ai < 



at < 



(^+2).(l-e) 



This implies that C 
for all e > 0, which yields 



k — e+2 2+log2fc- 

For a^, the proof is similar. The difference is that we now have three teams in every match, 
and two of the teams win. Let k = 3^. In the first round, {1} plays {2} and {3}, {4} plays {5} 
and {6}, and so on. In the second round, {1,2,3} plays {4,5,6} and {7,8,9}, and so on. We 
have only two final matches were {1, . . . , k/3} plays {k/3 + 1, . . . , 2A;/3} and {2A;/3 + 1, . . . , k}. 
Overall, we have £+1 matches. See Figure 1(b) for an example. Similar to Lemma 3.2, there 
exist now three (i — l)-underdogs, and they are in two different teams in the final matches. 
All three of them must win both finals in order to become a winner, which is impossible. This 
implies al<j^ = j^^- 

The lower bound states that every cycle cover without heavy-weight edges allows for a 
good decomposition. The proof uses the probabilistic method and Hoeffding's inequality [21, 
Theorem 2]. 

Lemma 3.3 (Hoeffding's inequality). Let Xi, . . . ,Xn be independent random variables, where 



Xj assumes values in the interval [aj,bj]. Let X = 



'{X < K{X) -t) < exp 



21^ 



Then 
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We consider the directed case first. Let A = 2 In A; + 4, and let C be an arbitrary directed 
cycle cover with it;(e) < 'w{C)/A for all e £ C. By scaling the edge weights, we make sure that 
Wi{C) = A and Wi{e) £ [0, 1] for all i £ [k]. 

Let ci,... ,Cm be the cycles of m. Consider any cycle Cj of C. We choose one edge of cj for 
removal uniformly at random. By doing this for j G [m], we obtain a decomposition P of C. 
Fix any objective i. Let Xj = ^eec nP^*(^) random variable of the contribution of cj 

to the weight Wi{P). Since Wi^e) £ [0, 1] for all e £ C, there exists aj,bj £ M such that Xj 
assumes only values in [uj, bj] and < bj — aj < 1. Let X = Yl^i -^j = Wi{P) be the random 
variable of the weight of P with respect to objective i. Let > 2 be the length of cycle Cj. A 
fixed edge of Cj is removed with a probability of Thus, 

E(X) = ^ E(X,) = E 1 - • "'^(^^•) ^ E 2 • ^^(^^-^ = 2 ■ = 2 • 

If we can show that P(X < 1) < then, by a union bound, P(3i £ [k] : < 1) < 1, 

which would imply the existence of a decomposition P with w{P) > 1 = w{C)/A and prove 
the theorem for directed graphs. Since < 6j — < 1, we have 

m mm m 

Wi{C) = A = Y^ Wi{cj) >^bj>Y^ bj - aj > ^{bj - ajf. 
j=i j=i j=i j=i 

We plug E{X) > ^ and t = ^ — 1 and J2jLiii>j ~ o^j)^ ^ ^ iiito Hoeffding's inequality and 
obtain 

/A 2 

exp(— In/c) = — . 

k 

The proof for undirected graphs goes along the same lines. The differences are that A = 
I • In A; + I and E{X) > ^. Thus, t > ^ - 1, and we obtain 

¥{X < 1) < exp - ^\ < exp -— + - = exp(-lnA:) 




which completes the proof of the theorem. □ 



3.2 Improved Decompositions 

The previous section showed that a-decompositions cannot, in general, yield more than a 

0(1/ log /c) fraction of the weight of a cycle cover. In order to obtain constant approximation 
ratios, independent of k, we have to generalize the concept of decompositions. 

Let C be a cycle cover, and let w = {wi, . . . ,Wk) be edge weights. We say that the pair 
{C,w) is 7-light for some 7 > 1 if w{e) < w{C)/^ for all e £ C. 

Theorem 3.4. Let e be arbitrary with < e < 1/2, and let k > 2 be arbitrary. Let C be a 

cycle cover, and let w = {wi, . . . ,Wk) be edge weights such that {C,w) is ^^^-light. 

If C is directed, then there exists a collection P ^ C of paths with w{P) > — e) • w{C). 
If C is undirected, then there exists a collection P C. C of paths with w{P) > (| — e) • w{C). 



8 



Proof. The proof uses Hoeffding's inequality. Wc start by considering the directed case. Let 
C be a directed cycle cover with edge weights w such that {C, w) is ^^"^ -light. With the same 
notation as in the proof of the lower bound of Theorem 3.1, we obtain that the probability of 
the event Wi{P) < — e) • Wi{C) for some fixed i is at most 

for k >2. The proof for undirected cycle covers is identical and therefore omitted. □ 

In fact, the proof above shows that if (C.w) is + 5)-light for any 5 > 0, then C admits 
a — e)- or (| — e) -decomposition: For any 6 > 0, the probability that Wi{P) is too small 
is strictly less than 1/k. We have used the (slightly worse) bound of to allow for an 
efficient randomized algorithm for finding decompositions. For the subsequent sections, we 
define rik,e = 

To conclude this section, we remark that no better decompositions are possible: First, even 
if we are content with an 0,(1/ log k) fraction of the weight, we need an 0(log /c)-light cycle 
cover. This follows from Theorem 3.1. 

Second, we cannot hope to get more than one half of the weight in case of directed cycle cov- 
ers: Consider a cycle cover that consists solely of cycles of length two and where all edges have 
equal weight. Then we lose one half of the weight no matter how we decompose. Analogously, 
an undirected cycle cover that consists solely of cycles of length three and where all edges have 
equal weight shows that we cannot hope to get more than a 2/3 fraction of the weight. To 
achieve better ratios than | and 5, techniques other than decomposition are needed. 

3.3 Finding Decompositions 

We know that decompositions exist due to Theorem 3.4. But, in order to use them in approxi- 
mation algorithms, we have to find them efficiently. In the remainder of this section, we devise 
two algorithms that do this job, a deterministic one and a faster randomized algorithm. 

The randomized algorithm is immediately obtained by exploiting Theorem 3.4: Assume that 
we have a cycle cover C with edge weights w such that {C,w) is l/r/fe^^-light. We randomly 
select one edge of every cycle of C for removal and put all remaining edges into P. The 
probability that P is not a(^ — e)-or (| — e) -decomposition (depending on whether C is 
directed or undirected) is bounded from above by 1/fc < 1/2. Thus, we obtain a decomposition 
with constant probability. By iterating this until a feasible decomposition has been found, we 
obtain a Las Vegas algorithm with expected linear running-time. 

Theorem 3.5. For every e > 0, every k >2, and every 1/rjk^^-light pair {C,w), where C is 
a directed cycle cover, a (5 — e) -decomposition can he computed by a randomized algorithm 

whose running-tim,e is expected to be linear in the input size. 

If C is an undirected cycle cover, a (| — decomposition can be computed analogously. 

For the deterministic algorithm, we assume again that wc have a l/r/^ g-light (C,w). We 
scale the weights such that Wi{C) = l/r]k,e for all i. Then w{e) < 1 for all e G C. The main idea 
is to find a kernel, i.e., to reduce an arbitrary instance to a new instance whose size depends 
only on k and e. 

First, wc normalize our cycle cover such that they consist solely of cycles of the shortest 
possible length. For directed cycle covers C, we can restrict ourselves to cycles of length two: 
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Any cycle c of length £ with edges ei, . . . , can be replaced by [i/2\ cycles (e2j_i, e2j) for 
j = 1, . . . , [^/2j. If i is odd, then we add an edge e^+i with w{e£^i) = and add the cycle 
(e£,e£_|_i). (Technically, edges consist of vertices, and we cannot simply reconnect them. What 
we mean is that we create new cycles of length two, whose edges have the same names and the 
same weights as the original cycles.) We do this for all cycles of length at least three and call 
the resulting cycle cover C . Now any decomposition P' of C yields a decomposition P of the 
original cycle cover C by removing the newly added edges e^_|_i if they are in P' . Furthermore, 
^i(e) < 1 for the new cycle cover C'. Analogously, undirected cycle covers can be normalized 
to consist solely of cycles of length three. 

Second, assume that we have two cycles c and c' in a normalized cycle cover with w{c) + 
w{c') < 1. Then we can combine c and c' to c: Let ei,e2 and e'i,e'2 be the edges of c and c', 
respectively. (For undirected cycles, we proceed analogously.) Then we can replace Cj and 
by Ci with w{ei) = w{ei) + w{e'j). The cycle cover plus edge weights thus obtained are still 
1/r/fc^g-light. We continue combining cycles greedily until no more combinations are possible. 

The resulting cycle cover contains at most 2k /rj^^^ cycles. Thus, an optimal decomposition 
can be found with a running-time that now only depends on k and e. The normalization can 
be implemented to run in linear time. 

Theorem 3.6. Fix £ > and k > 2. Then, for every l/rjk^^-light pair {C,w), where C is 
an directed cycle cover, a (5 — e) -decomposition can he computed by a deterministic algorithm 
whose running-time is polynomial in the input size. 

If C is an undirected cycle cover, a (| — e) -decomposition can be computed analogously. 

We call the procedure described above Decompose with parameters C, w, and e: (7 is a 
cycle cover (directed or undirected), w = {wi, . . . ,Wk) are k edge weights, and e > 0. Then 
Decompose(C, ■«;,£:) returns a (^ — e)- or (| — £)-decomposition P ^ C, provided that {C,w) 
is 77fe,£-light. 

4 Approximation Algorithms for MC-Max-TSP 

In this section, MaxCC-Approx denotes the randomized FPTAS for cycle covers. More 

precisely, let G be a graph (directed or undirected), w = (i«i, . . . ,tvi~) be edge weights, £ > 
and p G (0, 1]. Then MaxCC-Approx(G, uj, k, e,p) yields a (1 — £)-approximatc Pareto curve 
of cycle covers of G with weights w with a success probability of at least 1 — p. 

4.1 Multi-Criteria Max-ATSP 

Our goal is now either to use decomposition or to reduce the fc-criteria instance to a (A; — 
l)-critcria instance. To this aim, we put the cart before the horse: Instead of computing 
Hamiltonian cycles, we assume that they arc given. Then we show how to force an algorithm 
to find approximations to them. To obtain a 1/2 — £ approximate Pareto curve, we have to 
make sure that for every Hamiltonian cycle H, we have a Hamiltonian cycle H in our set with 
w{H)>{^^-e)-w{H). 

Fix £ with < £ < 2Hfc' ^ ^® Hamiltonian cycle, and let Pi = max{^i'j(e) | e G H} be 
the weight of the heaviest edge with respect to the ith. objective. Let /3 = P{H) = {Pi, . . . , Pk)- 
We will distinguish two cases. 
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In the first case, we assume that /? < (r/^ ^ — e'^) • w{H), i.e., H does not contain any heavy- 
weight edges. (Note that r/^^ — £^ > since rjk^^ = and e < 2iHjk-) modify our edge 
weights w to as follows: 




if wi^e) < P and 

if Wi{e) > Pi for some i. 



This means that wc set all edge weights exceeding (3 to 0. Since H does not contain any 
edges whose weight has been set to 0, we have w{H) = w^{H). Furthermore, for all subsets 
C of edges, we have w^{C) < w{C). The advantage of is that, if we compute a (1 — e) 
approximate Pareto curve of cycle covers with edge weights w^, we obtain a cycle cover to 
which we can apply decomposition to obtain a collection P of paths. This set P yields then a 
tour H that approximates H. This is stated in the following lemma. 

Lemma 4.1. Let e > be arbitrary. Let H be a directed Hamiltonian cycle with w{e) < 
{j]k,e ~ ■ w{H) for all e & H. Let P = P{H), and let be a {1 — e) approximate Pareto 
curve of cycle covers with respect to . 

Then contains a cycle cover C with w^{C) > (1 — s) ■ 'w{H) and w^{e) < rj^^^ ■ w^{C) for 
all e E C. This cycle cover C yields a decomposition P C. C with w{P) > — ^e) ■ w{H). 

Proof. Since the Hamiltonian cycle H is in particular a cycle cover, the set C*^ contains a cycle 
cover C with w^{C) > (1 — e) • w^{H) = (1 — e) • w{H). For every edge e G C and every i, we 

have wf{e) < {r]k,^ - e^) ■ w^{H) < ^^^^^ • u;f (C) < r/fc,^ • WiiC). The last inequality follows 
from rjk^e — < ^k,e ' (1 — £)) which is equivalent to > rjk^e ' £• Plugging in r]k,e = if^ yields 
> 2Hfe' which is valid and proves the preceding inequalities. 

The cycle cover C can be decomposed into a collection P C C of paths with w{P) > wf^{P) > 
(i - e) • w^{C). Thus, w{P) > {\- e) ■ {I - e) ■ w{H) > (5 - 2e) ■ w{H). □ 

In the second case, we assume that there exists an edge e = {u,v) e H and an i G [k] with 
Wi{e) > {r]k,e — • w{H). We put this edge into a set K of edges that we want to have in our 
cycle cover no matter what. Then we contract the edge e by removing all outgoing edges of u 
and all incoming edges of v and identifying u and v. (The resulting graph is still complete.) 
In this way, we obtain a slightly smaller Hamiltonian cycle H' = H \ {e}. Again, there might 
be an edge e' e H' and an i' G [k] with Wii{e') > {rjk^e - e^) ■ Wii{H'). (Since w{H') < w{H), 
there can be edges that are heavy with respect to H' that have not been heavy with respect 
to £[.) We put e' into K, contract e' and recurse. How long can this go on? There are two 
cases that can make an end: First, we might have obtained a Hamiltonian cycle H' that does 
not have any more heavy-weight edges, i.e., w{e) < (rj^^^ — " w{H') for all e G H'. In this 
case, we can apply Lemma 4.1 with decomposition. Second, we might have an i G [k] with 
Wi{K) > (i — e) • Wi{H), where H is our original Hamiltonian tour. Then we have collected 
enough weight with respect to the ith objective, and we can continue with only A; — 1 objectives. 
The next lemma gives an upper bound on the number of edges that have to be put into K. 

Lemma 4.2. After at most 

log(i+e) 



f{k,e)=k- 



log(l-%,£+£3) 
iterations, the procedure described above halts. 
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Proof. We prove that after at most f{k,e) iterations, there exists an i with Wi{K) > — e) • 
Wi{H). Let H = Hq, and obtained from H^^i by contracting one edge as described. Then, 
for every £, there exists an i such that Wi{H^) < (l — ??fc,£ + e^) • Wi{Hi-i). Thus, for £ = km, 
there exists an i with Wi{He) < (l — rik,e + e^)™ • Wi{H). If K contains km edges, then there 
exists an i with Wi{K) > Wi{H) • (l — (l — rjk^e + e^)"*). Now we have 

l-{l-Vk,e + eT>h-^ ^ {l-Vk,e + eT < i + ^ 

<^ m ■ log{l - i]k,s + < log(i + e) 

^ logd+e) 

4^ m > — — 

To complete the proof, we observe that both log(5 + e) and log(l — rjk^g + e^) are negative, 
the latter since rik,e > e^, which follows from rik,e = 2^ > the choice of £ < D 

To obtain an algorithm, we have to find /? and K. So far, we have assumed that we already 
know the Hamiltonian cycles that we aim for. But there is only a polynomial number of 
possibilities for /? and K: For all f3, we can assume that for all i there is an edge with 
Wi{e) = Pi- Thus, for every i there are at most 0{v?) choices for hence at most 0{'n?^) in 
total. The cardinality of K is bounded in terms of f{k,e) as shown in the lemma above. For 
fixed k and e, there is only a polynomial number of subsets of cardinality at most f{k,e). In 
fact, we can restrict ourselves to the subsets K that are path covers: A path cover is a subset 
K of edges such that K does not contain cycles and both the indegree and outdegree of every 
vertex is at most one. 

Overall, we obtain MaxATSP-Approx (Algorithm 1) and the following theorem. 

Theorem 4.3. For every k > 1, £ > 0, MaxATSP-Approx is a randomized \ — £ approxi- 
mation for k-criteria Max-ATSP whose running-time for a success probability of at least 1 — p 
is polynomial in the input size and log(l/p). 

Proof. We have to estimate three things: approximation ratio, running-time, and success prob- 
ability. All proofs arc by induction on k. For k = 1, the theorem holds since there is a deter- 
ministic, polynomial-time 2/3 approximation for mono-criterion Max-ATSP. In the following, 
we assume that the theorem is correct for A; — 1, all fixed e > 0, and all p> 0, and we prove it 
for k, fixed £ > 0, and all p > 0. 

Let us start by estimating the approximation ratio. For this purpose, we assume that all 
randomized computations are successful. Let H be an arbitrary Hamiltonian cycle. For a 
subset K C H, let Hk be H with all edges in K being contracted. Then, by Lemma 4.2, there 
exists a (possibly empty) set K C H of edges of cardinality at most f{k,e/2) with one of the 
two following properties: 

1. There exists an i with Wi{K) > (^ — |) • Wi{H). 

2. For all e G Hk, we have w{e) < (%,e/2 — (f)^) " '^{Hk)- 

In the first case, there exists an H' € "P-j^p with Wj{H') > (^ — ^)-Wj{HK) for all j G [fc] 
H' combined with K yields a Hamiltonian cycle H that satisfies w{H) > (5 — §) • w{H): First, 
we have 

WiiH)>w,iK)> (l-^) -wiH). 
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Vtsp ^ MAxATSP-AppROx(G,it;,A;,e,p) 
input: directed complete graph G = {V,E), k > 1, edge weig 

output: (5 — e) approximate Pareto curve "Ptsp for fe-C-Max-ATSP with a success probability 



of at least 1 — p 

1: if A; = 1 then 

2: compute a 2/3 approximation T-'tsp 

3: else 

4: for all subsets K C E with < f{k,e/2) such that K is a path cover do 

5: contract all edges of K to obtain Gk 

6: for all bounds (3 of {Gk,w) do 

7: Ck,I3 ^ MAxCC-ApPROx{Gk,w'^, k, |, ^„2fc+2Kfe,e/2) ) 

8: for all C G with w;^(e) < 77^,^/2 • i(;'^(C) for all e G C do 

9: P <— Decompose(C, w^^ , e/2) 

10: add edges to U P to obtain a Hamiltonian cycle H 

11: add H to Ptsp 

12: for all i <— 1 to A; do 

13: remove the ith objective from w to obtain w' 

14: r^^p ^ MaxATSP-Approx(Gk, A; - 1, f , ^^2k+ink.s/2) ) 

15: for all H' G V^^p do 

16: H ^KUH' 

17: add i7 to Ptsp 

Algorithm 1: Approximation algorithm for MC-Max-ATSP. 



Second, for j ^ i, we have 

wj{H) = wj{K) + Wj{H') > Wj{K) + - 1) • WjiHK) = " |) " ^ji^)- 

In the second case, let (3i = niax{wi{e) \ e G Hr} < {jlk,e/2 — (f)^) ' Wi{HK)- Then, 
according to Lemma 4.1, Ck,(} contains a cycle cover C with w{C) > (1 — |) • w{Hk) and 
u;^(e) < 77(t,£/2 • w{Hk)- Thus, C can be decomposed into a collection P of paths with w{P) > 

— |) • w{C) according to Lemma 4.1. Together with K, this yields a Hamiltonian cycle H 
with 

wiH) > w{P) + w{K) > - 1) • ^(C) + wiK) > Q - £ j . (1 - £) . ^(i/^) + wiK) 

= " I + t) ■ ^^^^^ + ^^^^ - " ^) ■ "^^^^^ 

Now we analyze the success probability and the running-time simultaneously. By the induc- 
tion hypothesis, the running-time of MaxATSP-Approx is polynomial in the input size and 
in log(l/p) for a success probability of at least 1 —p for A; — 1 objectives. 

Lines 4 and 6 yield at most Tn?^^'^^^^'^/'^^ iterations. Then we compute that many Pareto 
curves of cycle covers in line 7, each of polynomial cardinality. Some of these cycle covers are 
then decomposed to yield Hamiltonian cycles for "Ptsp- This can be implemented to run in 
polynomial time. 
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Lines 4 plus 12 yield k ■ n^^^^'^l'^) < 'n?k~^'^f^^'^/'^'> iterations. Each iteration requires a call of 
MaxATSP-Approx with k — 1 criteria, which is also polynomial. 

We have at most 2v?^'^'^^^^'^/'^^ calls of a randomized procedure, each called with an error 
bound of 2^2fc+2/(fe,e/2) ■ By a union bound, this gives an overall error probability of at most p. □ 

4.2 Multi-Criteria Max-STSP 

The {\- e) approximation for MC-Max-ATSP works of course also for MC-Max-STSP since 
the latter is a special case of the former. Our goal, however, is a ratio of (| — e). As a 
first attempt, one might just replace the (| — e)-decompositions by (| — e)-decompositions. 
Unfortunately, this is not sufficient since contracting the heavy-weight edges in undirected 
graphs is not as easy as it is for directed graphs: First, both statements "remove all incoming" 
and "remove all outgoing" edges are not well-defined in an undirected graph. Second, if we 
just consider all edges of one vertex as the incoming edges and all edges of the other vertex as 
the outgoing edges, we obtain a directed graph, which allows only for a ratio oi\ — e. 

To circumvent these problems, we do not contract edges e = {u, v}. Instead, we set the 
weight of all edges incident to m or u to 0. This allows us to add the edge e to any Hamiltonian 
cycle H' without decreasing the weight: We remove all edges incident to tt or u from H\ and 
then we add e. The result is a collection of paths. Then we add edges to connect these paths 
to a Hamiltonian cycle. The only edges that we have removed are edges incident to u or v, 
which have weight anyway. 

However, by setting the weight of edges adjacent to tt or f to 0, we might destroy a lot of 
weight with respect to some objective. To circumvent this problem as well, we have to consider 
larger neighborhoods of the edges va. K. In this way, we can add our heavy-weight edge (plus 
some more edges of its neighborhood) to the Hamiltonian cycle without losing too much weight 
from removing other edges. Lemma 4.4 justifies this. 

Before going into the details, let us fix some notation. Let H be an arbitrary Hamiltonian 
cycle. Let eo, ei, . . . , e^-i be the edges of in the order in which they appear in H (cq is chosen 
arbitrarily). Let ej = {vj,Vj-^i}, where arithmetic of the indices here and in the following is 
modulo n. 

Lemma 4.4. Let H be a Hamiltonian cycle as described above, let w = {wi, . . . ,Wk) be edge 
weights, and let ei, . . . ,ei be any i distinct edges of H. Then there exists a j & [£] such that 

w{ej)<'^.w{H). 

Proof. Suppose otherwise and assume without loss of generality that Wi{H) > for all i. We 
scale the weights such that Wi{H) = 1 for all i. Then for all j there is an ij with Wi.{ej) > 
J ■ Wi.{H) = J. Thus, Yl^j=iJ2i=iWi{ej) > X]j=i f ' Wi.{H) = k. But, since ah edges are 
distinct, we also have X]j=i Yli=i ^ Yl\=i = k, which is a contradiction. □ 

Our aim is now again to force the algorithm to find a Hamiltonian cycle H with w{H) > 
(I — e) • w{H). Like for the asymmetric TSP, we distinguish two cases. The first case is that 
H consists solely of light-weight edges, i.e., w{e) < (r)k,e — ■ w{H) for all e e H. We define 
P and in the same way as in the previous section. The undirected variant of Lemma 4.1 is 
similar to its directed counterpart. 
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Lemma 4.5. Let e > be arbitrary. Let H be an undirected Hamiltonian cycle with w{e) < 
{Vk,£ — • ui{H) for all e G iJ. Let (3 = (3{H), and let be a {1 — e) approximate Pareto 
curve of cycle covers with respect to w'-^ . 

Then C^^ contains a cycle cover C with w^{C) > (1 — e) ■ w{H) and w^{e) < rj^^^ ■ w'^{C) for 
all e G C. This cycle cover C yields a decomposition P C. C with w{P) > (| — 2£) ■ w{H). 

Proof. The proof is goes along the same lines as the proof of Lemma 4.1. We just have to 
replace 1/2 by 2/3. The details are therefore omitted. □ 

The second case is that H contains a heavy-weight edge. Without loss of generality, let 
this edge be cq. Then we put cq into our set K. However, we cannot contract edges since 
we have an undirected graph. Instead, we set the weight of all edges incident to vq and vi 
to 0. But in this way, we lose the weight of ei and e_i = e„_i. In order to maintain the 
approximation ratio, we have to avoid that we lose too much weight. Therefore, we consider 
paths that include cq- If we set the weight of all edges incident to vj, . . . , Vj' with j < < f 
to 0, we lose only the weight of ej-i and Cj'. To keep track of things, we also put ej, . . . , ej/_i 
into K. Furthermore, we put the two edges ej-i,ej', whose weight might get lost, into T. 
By Lemma 4.4, we can make sure that both ej-i and eji are not too heavy. Finally, we put 
Vj,. . . ,Vj' into the set L, which is the set of vertices whose edge weights have been set to 0. 
We denote the corresponding edge weights by w^, i.e., w^{e) = if e n L 7^ 0. 

Given any Hamiltonian cycle H\ we can now remove all edges of weight 0, add the edges 
Cj, . . . ,ej/_i, and connect the collection of paths to obtain a new Hamiltonian cycle H with 
w^{H) = w^{H'). The only edges that we cannot force to be in H are ej_i and ej/. In order 
to maintain a good approximation ratio we have to make sure that both arc light with respect 
to all objectives. This is were Lemma 4.4 comes into play: We choose j and / such that 
w{ej-i) + w{eji) < V ■ w{H) for some small u > Q that will depend on e and k and that we 
will specify later on. 

As we did for directed graphs, we recurse until cither Wi{K) > (| — e) • Wi{H) for some i 
or we have w^{e) < [rjk^e — ■ w^{H) so that we can apply Lemma 4.5 to and H. The 
following lemma is the counterpart for undirected graphs of Lemma 4.2. 



Lemma 4.6. Let 



g{k,e) = k ■ 



log(g+g) 



log(l - r]k,e + 

and V G (O, qt^^-^] • Then the procedure described above halts after at most g{k, e) iterations. 

Proof. Let w^^^ = w, and let w^^^ be the edge weights after the zth iterations. Similar to the 
proof of Lemma 4.2, wc show that after at most z < g{k,£) iterations, there exists an i with 

For every z, there exists an i with w'f\H) < (l — ?7a;,£ + ■ w^^~^\h). Thus, for z = km, 
there is an i with w^'\h) < (1 - r]k,e + e^)"" • Wi{H). 

Unfortunately and different to the proof of Lemma 4.2, we do not have w{K) = w{H) — 
w^^' (H) since some edges ended up in T. Instead, we only have It; (K) = w{H)-w('^\H)-w{T). 
We use the inequality w{T) < vzw{H) to bound the weight of T. Thus, 



w 



{K) > (1 - (1 - r}k,e + e^)"" - ^km) ■ w{H). 
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Vtsp ^ MAxSTSP-AppROx(G,it;,A;,e,p) 
input: undirected complete graph G = {V, E), k > 1, edge weig hts u; : S ^ N^ e > 
output: (I — e) approximate Pareto curve Vtsp for fe-C-Max-ATSP with a success probability 



of at least 1 — p 

1: if A; = 1 then 

2: compute a 61/81 approximation Vtsp 

3: else 

4: for all subsets K C E with \K\ < h{k,£/3) such that K is a path cover do 

5: let L be the set of vertices incident to K 

6: obtain from w by setting the weight of all edges incident to L to 

7: for all bounds /9 of {G,w^) do 

8: CL,f3 ^ MaxCC-ApPROx(G, W^^ , k, |, ^^2,+2Mfc,./3) ) 

9: for all C G Cl,^ with w^^{e) < ri^^^fs ■ w^l^{C) for all e G C do 

10: P ^ Decompose(C, tt;^^, e/3) 

11: remove edges of weight from P 

12: add edges to U P to obtain a Hamiltonian cycle H 

13: add H to T^TSP 

14: for alH ^ 1 to /c do 

15: remove the ith objective from to obtain w'^ 

16: ^ MAXATSP-ApPROX(G, w;'^, fc - 1, I, ,^,,+/,(fe,./3) ) 

17: for all H' G V^i^ do 

18: remove edges of weight from H' 

19: add edges to H' \J K to obtain a Hamiltonian cycle H 

20: add H to Ptsp 

Algorithm 2: Approximation algorithm for MC-Max-STSP. 



Now we set m = ^^^f' and show that the set K contains enough weight after km iterations: 
1 - (1 - r]k,e + - ykm > | - e ^ (1 

The implication holds since v < ^ Jj^ implies vkm < 1/6. □ 

We now have the key ingredients for our algorithm MaxSTSP-Approx (Algorithm 2). 
What remains to be done is to specify v more precisely and to analyze the approximation ratio 
as well as running-time and success probability of MaxSTSP-Approx. 

For z/, we assume that < < g:^^^-^, which is necessary for Lemma 4.6. In addition, the 
set T of edges that we cannot force to be included in our Hamiltonian cycles must not be too 
heavy. Otherwise, we cannot maintain our approximation ratio. We know that |r| < 2-g{k,£) 
by Lemma 4.6. For w{T) < e • w{H), we need w{e) < ^^^^ ■ w{H) for all e G T. By 

Lemma 4.4, this can be achieved for £ > |^HL£&£lj = h'{k,£). 

Let us summarize the above: In order to find the at most g{k, e) heavy edges, we need to 
examine subsets of edges of cardinality at most g{k,e). In order to keep T light, we have to 
consider at most h'{k,e) — 1 edges to either side of these edges. Thus, we have to examine 



3\m 



vkm < n + s 



(1 - Vk,e + £^ 



m > 



— log(l-?7fe ,+£3)" 
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subsets of edges of cardinality at most {2h'{k, e)) ■ g{k, e) = h{k, s). This number depends only 
on k and e, but not on the actual size of the input. 

Now wc arc prepared to prove the main theorem of Section 4.2. 

Theorem 4.7. For every k > 1, e > 0, MaxSTSP-Approx is a randomized | — e approxi- 
mation for k-criteria Max-STSP whose running-time for a success probability of at least 1 — p 
is polynomial in the input size and log(l/p). 

Proof. We will first estimate the approximation ratio given that all randomized computations 
are successful. After that, we will prove bounds for the error probability and the running-time. 
All proofs are by induction on k. For k = 1, the theorem holds since there is a deterministic, 
polynomial-time 61/81 approximation for mono-criterion Max-STSP. In the following, we as- 
sume that the theorem is correct for A; — 1, all fixed £ > 0, and all p > 0, and we prove it for 
k, fixed e > 0, and all p > 0. 

Let H be an arbitrary Hamiltonian cycle. According to Lemma 4.6 and the considerations 
above, there exists sets K,T C H such that \K\ -\- \T\ < h{k,£/3), if is a path cover and 
consists of at most g{k,£/3) paths, and T is the set of edges that connect K to the rest of H. 
Let L C V he the set of vertices incident to edges in K, and let be obtained from w by 
setting all edges incident to L to (this includes in particular all edges in K, which have both 
endpoints in L, and K, which have one endpoint in L). 

One of the following two properties holds: 

1. There exists an i with Wi{K) > (| — |) • Wi{H). 

2. For all eeH,we have w^{e) < {r]k,e/2, - (|)^) • w^{H). 

In the first case, there exists an H' G P^SP ^^i^^ w^{H') > (| - 1) -w^iH) for all j G [k]\{i}. 
We remove all edges from H' that have weight with respect to , which does not affect 
w^{H'). We obtain a Hamiltonian cycle H from K U H' (which is a collection of paths) by 
adding appropriate edges. Then we have 

Wi{H) > Wi{K) > - . Wi{H) 

and, for j / i, 

= (I - 1) ■ Mm - .,(r)) >(?-£).(!-£). ^,im >{l-e)- ^m- 

In the second case, let Pi = m.&x.{wf [e) \ e G H} < {rjk^^/^ — (e/3)^) • w^{H). Then, 
according to Lemma 4.5 (applied to and H), Cl,i3 contains a cycle cover C with w^^{C) > 
(l — I) • w^{H) and w^^{e) < g/3 • w^^{H). This cycle cover C can then be decomposed 
into a collection P of paths with w^{P) > (| — ^) • w^{H). We can assume that P does 
not contain any edges of weight 0. Together with K and some appropriate edges, this yields a 
Hamiltonian cycle H with 

wiH) > w{K) + w\P) > w{K) + _ 1^ . u,^{H) > - 1^ . {nj{H) - w{T)) 

.(M)■(-l)■»(^)KM)•'"<^)K^'-)■'"*^'• 
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What remains to be analyzed is the success probabihty and the running-time. By the 
induction hypothesis, the running-time of MaxSTSP-Approx is polynomial in the input size 
and log(l/p) for a success probability of at least 1 — p for — 1 objectives. 

Lines 4 and 7 yield at most ri'^k+2h{k,e/3) iterations. Then we compute that many approximate 
Pareto curves of cycle covers in line 8, each of polynomial cardinality. Some of these cycle covers 
are then decomposed to yield Hamiltonian cycles for T'tsp- This can be implemented to run 
in polynomial time. 

Lines 4 plus 14 yield at most kn^^^^'^^^^ < n'^k+2h{k,e/3) iterations. Each iteration requires a 
call of MaxATSP-Approx with k — 1 criteria, which is also polynomial. 

We have at most 2in?^^'^^''^^'^^'^^ calls of a randomized procedure, each called with an error 
bound of 2^2fc+2fe(fc,e/3) ■ By a union bound, this gives an overall error probability of at most p. □ 

5 Deterministic Approximations for 2-C-Max-STSP 

The algorithms presented in the previous section are randomized due to the computation of 
approximate Pareto curves of cycles covers. So are all approximation algorithms for MC- 
Min-TSP that we are aware of with the exception of a simple (2 -|- e) approximation for 
MC-Min-STSP [25]. 

As a first step towards deterministic approximation algorithms for MC-Max-TSP, we present 
a deterministic 61/243 ~ 0.251 approximation for 2-C-Max-STSP. The key insight for the 
results of this section is the following lemma. 

Lemma 5.1. Let M he a matching, let H he a collection of paths or a Hamiltonian cycle, and 
let w be edge weights. Then there exists a subset P C. H such that 

(i) PUM is a collection of paths or a Hamiltonian cycle (we call P in this case an M-feasihle 
set) and 

(ii) w{P) > w{H)/3. 

Proof. We prove the lemma by induction on \M\. For \M\ = 0, we can choose P = H and the 
lemma follows. 

Now let \M\ = i > and assume that the lemma holds for all smaller sets M and all H. We 
prove that it also holds for \M\ = i hy induction on \H\. The base case is -?/ = 0, for which 
the lemma obviously holds by setting P = H. Now assume that the lemma holds for \M\ = £ 
and \H\ < m and that we have \H\ = m. 

We distinguish two cases. The first case is that MCiH ^ 0. Then we set P = MDH and H' = 
H\M. By the induction hypothesis, there exists a P' H' such that w{P') > w{H')/3 and P' 
is an M-feasible set. Since P C M, we have P'uM = {P'UP)U M. Thus, also P = P' U P is 
an M-feasible set. Observing that w{P) = w{P') + w{P) > w{Hr\M) + w{H\M)/3 > w{H)/3 
completes the first case. 

The second case is that M D H = fj). Let e = argmax{t/;(e) | e G H} be a heaviest edge of H, 
and let /i, /2 G H he the two edges of H that are incident to e. Let H' = H \ {e, /i, /2}. (It 
can happen that /i or /2 do not exist, namely if H is not a Hamiltonian cycle but a collection 
of paths. But this is fine.) 

Let us first treat the case that e is incident to two edges zi,Z2 G M of the matching. (The 
matching M is not necessarily perfect.) Then we contract zi and Z2 to a single edge z that 
connects the two endpoints of zi and Z2 that are not incident to e and remove the two vertices 
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(a) (b) 

Figure 2: Contraction for the proof of Lemma 5.1: We keep e and remove /i,/2- Then we 
can contract zi and Z2 to z. 

incident to e (see Figure 2). Let M' = (M \ {zi,zi}) U {z}. Since e, /i, /2 are removed, iJ' and 
M' are a vahd instance for the lemma, i.e., M' is a matching and H' is a collection of paths 
{H' cannot be a Hamiltonian cycle). Wc can apply the induction hypothesis since |M'| < |M|. 

In this way, wc obtain an M'-feasible set P' C H' with w{P') > w{H')/3. Set P = P' U {e}. 
Since w{e) > w({e, /i, /2})/3, we have w{P) > w{H)/3. It remains to be proved that P is 
M-feasible: Every vertex in M U P has a degree of at most 2 by construction. Thus, the 
only possibility for P not to be M-feasible is that it contains a non-Hamiltonian cycle. By 
the induction hypothesis, M' U P' does not contain such a cycle. Thus, also M U P does not 
contain a cycle, since M U P is obtained from M' U P' just by replacing z by zi, e, and Z2- 

What remains to be considered is the case the e is not incident to two edges zi,Z2 G M. 
Then we consider the shortest path in ei, . . . , € of edges in H that includes e such that 
ei and arc incident to any edges zi,Z2 G M. The reasoning above holds in exactly the 
same way if replace e by the path ei, . . . , Cg, and we put ei, . . . , into P. If no such path 
exists, then either M = 0, which we have already dealt with, or the path containing e ends 
somewhere at a vertex of degree 1 in H U M. In the latter case, we can simply put the whole 
path into P. □ 

Lemma 5.1 yields tight bounds for the existence of approximate Pareto curves with only a 
single element. This is the purpose of the following theorem. 

Theorem 5.2. (a) For every undirected complete graph G with edge weights wi and W2, there 
exists a Hamiltonian cycle H such that is a 1/3 approximate Pareto curve for 2-C- 
Max-STSP. 

(h) Part (a) is tight: There exists a graph G with edge weights wi and W2 such that, for all 
e > 0, no single Hamiltonian tour of G is a (1/3 + e) approximate Pareto curve. 

Before embarking on the proof of the theorem, let us make some remarks. The question 
how well a single element can approximate a whole Pareto curve has also been addressed by 
Ehrgott [11]: He proved that for MC-Min-STSP, there is always a single Hamiltonian cycle H 
such that the norm of w{H) is at most twice as large as the norm of w{H) for any Hamiltonian 
cycle H, i.e., ||tt;(ff)|| < 2 • ||tt;(i7)||. This, however, does not imply that {H} is already an 
approximate Pareto curve. Instead, single-element approximate Pareto curves exist for no 
other variant of multi-criteria TSP than 2-C-Max-STSP: For A;-C-Max-STSP for A; > 3, we can 
consider a vertex incident to three edges of weight (1,0,0), (0,1,0), and (0,0,1), respectively. 
All other edges of the graph have weight 0. Then no single Hamiltonian cycle can have positive 
weight with respect to all three objectives simultaneously. Similarly, no such result is possible 
for fc-C-Max-ATSP and for /c-C-Min-TSP for any k>2. 
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Figure 3: The graph for Theorem 5.2(b). Sohd edges are of weight (1,0), dashed edges are 
of weight (0, 1), all other edges (dotted or not drawn) have weight (0,0). 

Proof. Let Hi and H2 Hamiltonian cycles of G such that Hi maximizes wi and H2 maximizes 
W2- Then there exists a matching M C Hi with w{M) > w{Hi)/3. (We can actually get 
w{Hi)/2 if G has an even number of vertices and ■ w{Hi) if the number n of G's vertices 
is odd. This, however, docs not improve the result.) We apply Lemma 5.1 with H = H2 and 
obtain an M-feasible set P C H2. From M and P, we obtain a Hamiltonian cycle H' D MUP: 
Either M U P is already a Hamiltonian cycle, then nothing has to be done. Or M U P is a 
collection of paths. Then we add appropriate edges to obtain H'. We claim that {H'} is a 1/3 
approximate Pareto curve: Let H be any Hamiltonian tour. Then 

wi{H') > wi{M) > wi{Hi)/3 > wi{H)/3 

and 

W2{H') > W2{P) > W2{H2)/3 > W2{H)/3. 

To finish the proof, let us show that the 1/3 bound is tight, i.e., in general, no single solution 
can be better than a 1/3 approximate Pareto curve. The proof is by contradiction. Consider 
the graph shown in Figure 3: The solid edges have weight (1, 0), the dashed edges have weight 
(0, 1). All other edges (two are dotted, the other edges are not shown) are of weight (0, 0). 

The solid edges plus the dotted edges form a Hamiltonian cycle of weight (3, 0). The dashed 
edges plus the dotted edges form a Hamiltonian cycle of weight (0, 3). Assume that there exists 
a Hamiltonian cycle H such that {H} is a 1/3 + e approximate Pareto curve for our instance. 
Then H must contain at least two solid edges (otherwise wi{H) < 1) and at least two dashed 
edges (otherwise W2{H) < 1). In order to include the top vertex, we also need two edges of 
weight (0,0). Thus, in total at least six edges are needed, but any Hamiltonian cycle on a 
five-vertex graph contains only five vertices - contradiction. □ 

Lemma 5.1 and Theorem 5.2 are constructive in the sense that, given a Hamiltonian cycle 
H2 that maximizes W2, the tour H can be computed in polynomial time. A matching M with 
wi{M) > 'Wi{Hi)/3 can be computed in cubic time. However, since we cannot compute an 
optimal H2 efhciently, the results cannot be exploited directly to get an algorithm. Instead, 
we use an approximation algorithm for finding a Hamiltonian tour with as much weight with 
respect to W2 as possible. Using the 61/81 approximation algorithm for Max-STSP [10], we 
obtain Algorithm 3 (which in particular is an algorithmic version of Lemma 5.1) and the 
following theorem. 

Theorem 5.3. BiMaxSTSP-Approx is a deterministic 61/243 approximation algorithm 
with running- time 0{n^) for 2-C-Max-STSP. 

Proof. The running-time is dominated by the running-time of the 61/81 approximation for 
Max-STSP by Chen et al. [10] and the time for computing the matching, both of which is 
O(n^). The approximation ratio follows from |j • | = □ 
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Vtsp ^ BiMAxSTSP-AppROx(G,u;i,ti;2) 


input: undirected complete graph G = {V, E), edge weights wi,W2 '■ E N'' 


output: a 61/243 approximate Pareto curve H 


1 


compute a maximum-weight matching M with respect to wi 


2 


compute a 61/81 approximate Hamiltonian cycle H2 with respect to W2 


3 


P ^ H2r\M 


4 


M' ^ M 


5 


Ho^ HoXP 


6 


while il2 7^ do 


7 


6 < — argmax{'U?2 (^0 | ^ -^^2} 


8 


pxtPTiH p to fi natVi Pi p^ ffn snrli tliat, nnl v Pi arifl arp inriHpnf, t,o pHp'ps 
, Z2 G M' or the path cannot be extended anymore 


9 


P^PU{ei,...,eg} 


10 


H2^ H2\ {ei, . . . ,eg} 


11 


if zi or Z2 exists then 


12 


let /i , /2 € -^2 be the two edges extending the path if they exist 


13 


^^2^i^2\{/l,/2} 


14 


if both zi and Z2 exist then 


15 


contract zi and Z2 to z 


16 


M' ^{M'\{zuZ2})V^{z} 


17 


let be a Hamiltonian cycle obtained from P U M 



Algorithm 3: Approximation algorithm for 2-C-Max-STSP. 



For the special case of metric 2-C-Max-STSP, i.e., the edge weights are restricted to fulfil 
the triangle inequality, we obtain the slightly better approximation ratio of 7/24 > 0.29 if we 
replace the 61/81 approximation with the 7/8 approximation for metric Max-STSP by Kowalik 

and Mucha [24]. 

Finally, we note that the bound of 1/3 of Theorem 5.2 holds only for small instances. If 
the weight of the optimum Hamiltonian cycles is distributed among many edges, then (| — 5) 
approximate Pareto curves exist for some (5 > 0. The following lemma states this more precisely. 

Lemma 5.4. Let G = {V^ E) he an undirected complete graph, and let wi and W2 be edge 
weights. Let Hi and H2 be two Hamiltonian cycles maximizing wi and W2, respectively. Let 
c = max| J^'^^^^ I i G {1, 2}, e G Hi\ be the maximum contribution of a single edge to w\{Hi) 
or W2{H2). 

Then there exists a Hamiltonian cycle H that is a single- element (| — approximate Pareto 
curve for 5 > ^ + ■ In 2. 

Proof. We prove that there exists a collection P of paths such that Wi{P) > (| — 5) • Wi{Hi) 

for z = 1,2. There exists matchings Mj C Hi with Wi{Mi) > ""'^^^^""^ for i G {1,2}. Consider 

the graph (V, Mi U M2). Any component of this graph is either a single edge that is common 
to Ml and M2. Then we put this edge into P. Or it is a cycle of even length, consisting in 
turn of edges of Mi and M2. Then we have to remove an edge to obtain a path. Let ci, . . . , q 
be these cycles. 

For every cycle, we choose randomly to remove either the lightest edge from Mi or the 
lightest edge of M2. Thus, the contribution of Cj to Wi is a random variable Xj that assumes 
values in an interval of at most c. We have lE(X]j=i Xj) > | • Wi{Mj) > ^Wi{Hi) — ^. Thus, 
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by Hoeffding's inequality, the probability that ^{Yl^j=i -^j) < (1 ~ ' "^iiHi) is at most 

where aj is Wi[cj) minus the weight of the hghter edge of Mi with respect to Wi and hj = Wi{cj). 
We have bj — aj < c and Yl^j=i = 1 as well as bj < 2a j , which yields 



1 ^ 



Thus, 



Altogether, 



2 



c 



5>- + ^-.ln2. 

Thus, by a union bound, the probability that wi{P) < | — (5 or W2{P) < | — (5 is smaller than 
1, which implies the existence of such a P. □ 

For instance, if c < 0.01, i.e., no edge contains more than 1% of the weight of an optimal 
Hamiltonian cycle, then a single-element 0.37 approximate Pareto curve exists. If c < 0.18, 
then a single-element 0.35 approximate Pareto curve exists. 



6 Approximation Algorithm for MC-Min-ATSP 

Now we turn to MC-Min-ATSP and MC-Min-7-ATSP, i.e., Hamiltonian cycles of minimum 
weight are sought in directed graphs. Algorithm 4 is an adaptation of the algorithm of Frieze 
et al. [16] to multi-criteria ATSP. Therefore, we briefly describe their algorithm: We compute 
a cycle cover of minimum weight. If this cycle cover is already a Hamiltonian cycle, then 
we are done. Otherwise, we choose an arbitrary vertex from every cycle. Then we proceed 
recursively on the subset of vertices thus chosen to obtain a Hamiltonian cycle that contains 
all these vertices. The cycle cover plus this Hamiltonian cycle form an Eulerian graph. We 
traverse the Eulerian cycle and take shortcuts whenever visiting vertices more than once. The 
approximation ratio achieved by this algorithm is log2 n for Min-ATSP [16] and 1/(1 — 7) for 
Min-7-ATSP [7]. 

To approximate MC-Min-ATSP, we use MinATSP-Approx (Algorithm 4), which pro- 
ceeds as follows: We compute an approximate Pareto curve of cycle covers. This is done 
by MinCC-Approx, where MinCC-Approx(G, it;, A;, e,p) computes a (1 + e) approximate 
Pareto curve of cycle covers of G with weights w with a success probability of at least 1 — p 
in time polynomial in the input size, 1/e, and log(l/p). (Of course, the aim is now to find 
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Vtsp ^ MinATSP-Approx(G, w, k, e) 
input: directed complete graph G = {V,E) with n = \V\, k > 1, edge weights w : E N*^, 
£ > 

output: (logn + e) approximate Pareto curve for A;-C-Min-ATSP or + e) approximate 
Pareto curve for /c-C-Min-7-ATSP with a probabihty of at least 1/2 

1: e' ^ eV log^ n; ^ 0; j ^ 1 

C ^ MinCC-Approx(G, w, k, e', 2Q\ogn ) t> Q is defined in Lemma 6.2 

Vo^{{C,w{C),y,-^)\CeC} 
while ^ do 

for all TT = (C, w', V, n') G Vj-i do 
if {V, C) is connected then 
^ J^U{{C',w'y,TT')} 

else 

select one vertex of every component of (F', C) to obtain V 
C ^ MiNCC-AppROx(G,i(;,A;,e', 2Q^) 
Vj ^ Vj U {((7, -u), F, tt) I (7 G C, tt) = u;' + 7^' • u;(C')} 
while there arc 7r',7r" G with the same e'-signature do 
remove one of them arbitrarily 

J ^ J + 1 
T'tsp ^ 

for all {C\w',V'y) G do 
H 

while vr' = {C\ w\ V\ tt") ^ ± do 

construct a Hamiltonian cycle H' on V" from H U C" by taking shortcuts such that 
all edges of H are removed 

TT <— TT 



H ^ H' 
Vtsp ^ Vtsp U {H} 



Algorithm 4: Approximation algorithm for MC-Min-ATSP and MC-Min-7-ATSP. 



cycle covers of minimum weight.) Then we iterate by computing approximate Pareto curves of 
cycle covers on vertex sets V for every cycle cover C in the previous set. The set V contains 
exactly one vertex of every cycle of C. Unfortunately, it can happen that we construct a super- 
polynomial number of solutions in this way. To cope with this, we remove some intermediate 
solutions if there arc other intermediate solutions whose weight is close by. We call this process 
sparsification. It is performed in lines 13 and 14 of Algorithm 4 and based on the following 
observation: Let e > 0, and consider H of weight w{H) G N*^. For every i G {1, . . . , fc}, there 
is a unique £i G N such that Wi{H) G [(1 + e)^' , (1 + £)^'+^) . We call the vector £ = (£1 , . . . , 4) 
the e-signature of H and of w{H). Since w{H) < 2^^^\ where N is the size of the instance, 
ii is bounded by a polynomial q{N, 1/e). There are at most different e-signatures, which is 
polynomial for fixed k. To get an approximate Pareto curve, we can restrict ourselves to have 
at most one solution with any specific e-signature. 

In lines 2 and 3 of Algorithm 4, an initial (1 + e') approximate Pareto curve of cycle covers 
is computed. In the loop in lines 4 to 15, the algorithm computes iteratively Pareto curves 
of cycle covers. The set Vj contains configurations vr = (C, w', F', tt'), where C is a cycle 
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Figure 4: Comparison of the approximation ratio subject to 7 for MC-Min-7-ATSP of Man- 
they and Ram's algorithm [25] and our algorithm. 

cover on V' , n' is the predecessor configuration, and w' is the weight of C plus the weight 
of its predecessor cycle covers, each weighted with an appropriate power of 7. (We define 
the e'-signature of tt = {C ,w' ,V' ,Tr') to be the e'-signature of w'.) The reason for these 
weights will become clear in the analysis of the approximation ratio of the algorithm. If, in 
the course of this computation, we obtain Hamiltonian cycles, these are put into (line 8). 
In lines 13 and 14, the sparsification takes place. Finally, in lines 16 to 23, Hamiltonian cycles 
are constructed from the cycle covers computed. 

Let us now come to the analysis of the algorithm. Our goal is to prove the following 
result, which follows from Lemmas 6.2, 6.3, and 6.6 below. MinATSP-Approx is the first 
approximation algorithm for MC-Min-ATSP and for MC-Min-7-ATSP for 7 > ^/TJ3 « 0.58. 
Furthermore, for 7 > 0.55, it improves over the previously known algorithm, which works only 
for 1/2 < 7 < \/l73 0.58 (see Figure 4 for a comparison of the approximation ratios). 

Theorem 6.1. For every e > 0, Algorithm 4 is a randomized (logn + e) approximation for 

MC-Min-ATSP and a randomized +e) approximation for MC-Min-j-ATSP for 7 G [|, 1) 
with a success probability of at least 1/2. Its running-time is polynomial in the input size and 
1/e. 

Wc observe that for every j and every vr = (C ,w' ,V' ,7r') G Vj, we have \V'\ < For 
j = 0, this holds by construction (line 3). For j > and n' = {C" , V" , w" ,7r") G Vj-i, we have 
< by the induction hypothesis. Since every cycle involves at least two vertices, 

we have \V'\ < \V"\/2 < This yields also that Vj is empty for j > [lognj: Assume to 

the contrary that such a Vj contained a configuration {C ,V ,w,7r). Then \V\ < n/2Li°g"J ^ 2. 
Thus, \V\ < 1. Let tt = (C ,w' ,V' ,-k'), then this implies that {V',C') had been connected. 
Hence, we would enter line 8 rather than lines 10 to 12 - a contradiction. 

Let us now analyze the running-time. After that, we examine the approximation performance 
and finally the success probability. 

Lemma 6.2. The running-time of Algorithm 4 is polynomial in the input size and 1/e. 

Proof. Let N be the size of the instance at hand, and let Q = Q{N, 1/e') be a two-variable 
polynomial that bounds the number of different e'-signatures of solutions for instances of size 
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at most N. We abbreviate "polynomial in the input size and simply by "polynomial." 
This is equivalent to "polynomial in the input size and 1/e'" by the choice of e'. 

The approximate Pareto curves can be computed in polynomial time with a success proba- 
bility of at least 1 — {2Q\ogn)~^ by executing the randomized FPTAS \log{2Q\ogn)] times. 
Thus, all operations can be implemented to run in polynomial time provided that the cardi- 
nalities of all sets Vj are bounded from above by a polynomial Q for all j. Then, for each j, 
at most Q approximate Pareto curves of cycle covers are constructed in line 11, each one in 
polynomial time. 

For every e'-signature and every j, the set Vj contains at most one cycle cover with that spe- 
cific e'-signature. The lemma follows since the number Q of different e'-signatures is bounded 
by a polynomial. □ 

Let us now analyze the approximation ratio. First, we will assume that all randomized 
computations of (1 + e') approximate cycle covers are successful. After that we analyze the 
probability that one of them fails. 

Lemma 6.3. Assume that in all executions of line 11 of Algorithm ^ an (1 + e') approxi- 
mate Pareto curve of cycle covers is successfully computed. Then Algorithm 4 achieves an 
approximation ratio o/X^y°o""'~^ 7'' + ^ /^^^ MC-Min-j-ATSP for 7 G [5, 1]. 

If 7 < 1, then ZfT^'' 1' < J:T=o 7^' = t^. If 7 = 1, then E^"^"' 7^' = [log nj < log n. 
Thus, we obtain the approximation ratios claimed in Theorem 6.1. 

Proof Let r = Ei=r^~S^ and let H be any Hamiltonian cycle on V. We have to show 
that the set Ptsp of solutions computed by Algorithm 4 contains a Hamiltonian cycle H with 
w{H) <{r + £)- w{H). 

Lemma 6.4. Let tt = {C ,w' ,V' ,Tr') G T. Then, in lines 16 to 23, a Hamiltonian cycle H is 
constructed from it with w{H) < w' . 

Proof. Assume that tt^ = tt for some i and that ttj = {Cj,Vj,w^\Trj-i) for j G {0,... 
(We have 7r_i = _L, C = Q, V = Ve, u;^ = w', 7r£_i = vr', and Vq = V.) Let He = Q be 
the Hamiltonian cycle on V^. In lines 16 to 23, we iteratively construct Hamiltonian cycles Hj 
on Vj from Cj and i^j+i- Wc have w{Hj) < w(Cj) + 7 • w{Hj^i) by the 7-triangle inequality. 
Thus, w{H) = w{Ho) < EJ=o 7^ ' ^(^i) = ^' by definition of w'. □ 

What remains to be proved is that, for every Hamiltonian cycle H, there exists a tt = 
(C, w', V, it') in JF such that w' < {r + e) ■ w{H). 

Lemma 6.5. For every £, there exists a configuration n = {C ,vj' ,V' ,7r') G Vi U with 
w' <{l + s'Y+^-Y!j=ol' -^{H). 

Proof. The proof is by induction on For £ = 0, the lemma boils down to the existence of a 
configuration (C, u;', F', _L) G Vq with w{C') < (1 -|- e') ■ w{H). Such a C exists because in 
line 2, a (1 -|- e') approximate Pareto curve of cycle covers is computed. 

Now assume that the lemma holds for — li T contains a configuration (C, w' , V, tt') that 
satisfies the lemma for £ — 1, then we are done since (1 + s'Y ■ Ej=o 7"' ^ (1 + s'Y~^^ ■ Ej=o 7"'- 
Otherwise, Ve-i contains a configuration tt' = (C", V" , w" ,tt") with w" < (1 + e'Y • Ylj=o 7"' ' 
w{H). Let V be the set of vertices constructed from tt' in line 10, and let H' be H restricted 
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to V by taking shortcuts. By the triangle inequahty, we have w{H') < w{H). After hne 11, 
C contains a cycle cover C with w{C') < (1 + e') • w{H'). Let tt = {C ,w' ,V' ,11') with 
w' = w" + -wiC). Then 

w' < ^{1 + e'Y • ^ 7^' + 7^ • (1 + e') j • wiH). 

What remains to be analyzed is the sparsification in lines 13 to 14. If tt G Vj after sparsification, 
we are done since 

£—1 i 

(1 + e'Y ■ y + (1 + e') • 7^ < (1 + sT' • E 7^'- 

j=0 j=0 

Otherwise, Ve contains a tt = {C' ,w' ,V' ,Tt') with the same e'-signature as tt. Thus, 

i-i 



w 



' < (1 + s') ■[{1 + e'Y + e') ■ w{H) < (1 + • E " ^(^)' 

j=0 / j=0 



and TT fulfills the requirements of the lemma. □ 

Let r be defined as in the beginning of the proof of Lemma 6.3. By Lemmas 6.4 and 6.5, 
we obtain an approximation ratio of (1 + £')L'°snJ . j. since Vg, is empty for i > [lognj. What 
remains to be proved is that this is at most r + e: 

r • (1 + £')Li°g"J <r-(l + -ir-] ^ < r • exp f ^J-^ <r-(l + <r + e. 

V log'^n/ Vlog n/ V logny 

The first inequality follows from our choice of e'. The second inequality holds since (l + < 
exp(x). The third inequality holds because exp(x^) < 1 + x for a; G [0,0.7] (we assume 
e/logn < 0.7 without loss of generality.) The fourth inequality holds since r < logn. □ 

Lemma 6.6. The probability that, in a run of MinATSP-Approx (Algorithm 4), i'n the 
execution of every line 2 and in every execution of line 11 an e' approximate Pareto curve of 
cycle covers is successfully computed is at least 1/2. 

Proof. Lines 11 and 2 of Algorithm 4 is executed at most Q ■ logn times, where Q is an upper 
bound for the number of different e'-signatures of solutions of instances of size at most N . (Q 
is polynomial in N and l/e.) Each execution fails with a probability of at most l^g „ ■ Thus 

by a union bound, the probability that one of them fails is at most 2Q^°og'w = 1/2- D 

MinATSP-Approx uses randomization only for MinCC-Approx. Thus, Lemma 6.6 im- 
mediately yields that our algorithm has a success probability of at least 1/2 in compliance 
with Theorem 6.1 and the definition of a randomized approximation algorithm. The success 
probability can be amplified to 1 — p for arbitrarily small p > by running the algorithm 
0(log(l/p)) times. 
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7 Conclusions 



We have presented approximation algorithms for almost all variants of multi-criteria TSP. The 
approximation ratios of our algorithms are independent of the number k of criteria and come 
close to the currently best ratios for TSP with a single objective. Furthermore, all of our 
algorithms work for any number of criteria. 

In particular, we have presented a factor ^ — £ approximation algorithm for MC-Max-ATSP 
and a factor |— e approximation algorithm for MC-Max-STSP. The algorithms arc randomized, 
and their running-time is polynomial in the input size for all fixed e > and k. In terms of 
approximation ratios, this is close to the best known ratios for mono-criterion TSP (| for Max- 
ATSP and |]- for Max-STSP). Our algorithms improve upon the previous algorithms for these 
problems that achieve ratios of l/(fc -|- 1) -|- £ and 1/k + e, and they give an affirmative answer 
to the question raised by Blaser et al. [6] whether there exist algorithms with approximation 
ratio $7(1/ log fc). 

Furthermore, we have presented a randomized log2 n + s approximation for MC-Min-ATSP. 
The approximation ratio is asymptotically equal to the best known ratio of | • log2 n for Min- 
ATSP. Finally, wc devised a deterministic 61/243 approximation for 2-C-Min-STSP with cubic 
running-time, and we proved that for 2-C-Min-STSP, there always exists a 1/3 approximate 
Pareto curve that consists of a single element. 

Most approximation algorithms for multi-criteria TSP use randomness since computing ap- 
proximate Pareto curves of cycle covers requires randomness. This raises the question of 
whether there are algorithms (particularly designed) for multi-criteria TSP that are faster, 
deterministic, and achieve better approximation ratios. 
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